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Введение. Верифицирована методика численного модели- 
рования поперечного обтекания пролетных строений мо- 
стов на основе нестационарного решения Рейнольдса для 
уравнений Навье — Стокса (ОКАМ$, Опуеа4у Кеупо!&$- 
ауегаге4 Ма\мег — зюКез). Данный двумерный подход 
используется в современных методах и пакетах приклад- 
ных программ вычислительной гидроаэродинамики. Цели 
работы — отладка и экспериментальное обоснование ука- 
занной методики. Для реализации поставленной цели ис- 
пользована ранее разработанная авторами база данных по 
аэродинамическим характеристикам поперечных сечений 
пролетных строений балочных мостов типовых форм. 
Материалы и методы. Проведено численное моделирова- 
ние поперечного обтекания мостовых строений низкотур- 
булентными (гладкими) и турбулентными воздушными 
потоками в диапазоне практически интересных углов ата- 
ки. В итоге использовалась модель турбулентности Э5Т 
К — «. Методика предварительно отработана на тестовой 
задаче обтекания балок прямоугольного поперечного сече- 
ния. Расчеты проводились с помощью лицензионного про- 
граммного комплекса АМЗУЪ. 

Результаты исследования. Показано, каким образом угол 
атаки определяет силы (подъемную и лобового сопротив- 
ления) и момент поперечных сечений балочных мостов 
типовых форм. Полученные расчетные зависимости отно- 
сятся к пролетным строениям на стадиях монтажа (без 
плиты перекрытия и ограждений, без ограждений) и экс- 
плуатации в условиях модельных набегающих потоков — 
гладкого и турбулентного. Это позволяет очертить грани- 
цы для более взвешенных оценок аэродинамических ха- 
рактеристик балочных мостов в реальном ветровом потоке. 
Лобовое сопротивление сечению демонстрирует наилуч- 
шее согласование с опытными данными. Величина подъ- 
емной силы более чувствительна к наличию и протяжен- 
ности отрывных зон, поэтому ее расчетное определение 
менее точно. Наиболее проблемным для большинства 
конфигураций является воспроизведение влияния угла 
атаки на аэродинамический момент сечения. 
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Обсуждение и заключения. Сравнение расчетных и опыт- 
ных данных свидетельствует о применимости ОКАХ$- 
подхода к оперативному предсказанию аэродинамических 
характеристик однобалочных пролетных строений. Если 
же речь идет о многобалочных пролетных строениях, су- 
щественную роль играет аэродинамическая интерференция 
между отдельными балками. В этом случае вместо 
ОВАЮ-подхода следует применять более точные вихре- 
разрешающие методы. Полученные результаты могут быть 
использованы в исследованиях аэродинамики сооружений 
и в практике проектных организаций в сфере транспортно- 
го строительства. 
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стов; аэродинамические характеристики. 
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Введение. Известно, что обтекание ветром инженерных сооружений носит, как правило, нестационар- 
ный турбулентный характер, в потоке наблюдаются разномасштабные вихревые структуры [1, 2]. Вблизи пло- 
хообтекаемых тел (к которым, в частности, относятся пролетные строения мостов) возникают нестационарные 
отрывные зоны [3]. В связи с этим адекватное моделирование эффектов турбулентности становится важным 
требованием к методике вычислительного эксперимента. 

Для оценки аэродинамики сооружений используется моделирование крупных вихрей (Гагое Еа4у 
эпишаНоп, ГЕ) и гибридное описании потока (Реасвеа Едду Зипайоп, РЕЗ). Однако применение указанных 
методов осложнено их повышенной ресурсоемкостью, причины которой: 

— трехмерность задачи; 

— высокие требования к плотности расчетной сетки в пристеночной области и в области «фокуса» [4]; 

— ограничения на шаг интегрирования по времени; 

— относительно большая длительность временного «окна» для сбора нестационарной статистики в установив- 
шемся режиме. 

В то же время известно [5], что для цилиндрических призм, близких по форме к балочным мостовым 
пролетным строениям, двумерный подход воспроизводит основные свойства течения (первичная неустойчивая 
мода в следе за телом существенно двумерна). Так, в работе [6] выполнен РОЛ-анализ (Ргорег Огфогопа| Ое- 
сотроз1Яоп — истинное ортогональное разложение) течения вблизи призмы с относительной шириной сечения 
В/Н =5 (Н — высота сечения). В результате установлено, что 1-я и 2-я моды возмущений двумерны (постоян- 
ны по пролету) и соответствуют переносу вихрей вдоль поверхности. Трехмерные моды изменяются по пролету 


на характерной длине, которая не меньше ширины сечения Б. 


Для оперативного предсказания аэродинамических характеристик (АДХ) мостовых сооружений и со- 
провождения продувок в аэродинамической трубе авторы настоящей работы использовали нестационарное 
двумерное моделирование на основе осредненных по Рейнольдсу уравнений Навье — Стокса (ОКАХ-подход, 
Опзеаду Кеупо!4$-ауегагеа Мау1ег — 5юКез$). Его применимость к определению АДХ плохообтекаемых тел 


(неподвижных и колеблющихся) изучалась в ряде работ зарубежных авторов (см., например, [7]). 


Механика 


Материалы и методы. При постановке компьютерных экспериментов учитывались приведенные в [8— 


10] рекомендации. Расчеты проводились в программе Ншеп( комплекса АМЗУЪ. 
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Предварительно методика была отработана на тестовой задаче обтекания балок прямоугольного попе- 
речного сечения. В итоге для описания турбулентности потока была выбрана модель переноса сдвиговых 
напряжений (ЭВеаг 5гез5 Тгапзрой, э5Т) К — ®, подобраны параметры сетки и численного алгоритма. 

Протяженность прямоугольной расчетной области в продольном направлении — (30...40)Н, в попе- 
речном — (14... 20) Н. Передняя грань обтекаемого тела отстояла от входной границы на (8 ...12)Н. 

Использовались низкорейнольдсовые сетки (безразмерное расстояние первого узла до стенки у* < 
4...5), позволяющие рассчитывать отрыв и присоединение пограничного слоя. Ввиду сложности контуров об- 
текаемых тел строились многоблочные сетки. Внутренний, охватывающий тело, блок состоял из четырехуголь- 
ных элементов, плотность которых ближе к поверхности тела увеличивалась. Непосредственно у стенки созда- 
вался слой со структурированной ортогональной четырехугольной сеткой. Область следа покрывалась сеткой 
из квадратных ячеек размером не более Н/15...Н/10. К внешним границам размер ячеек возрастал до Н/ 
4...Н/З. На периметр сечения в зависимости от его формы приходилось порядка 10°...10? ячеек. Общее коли- 
чество ячеек варьировалось от 40—50 тысяч (для сечений простых форм) до 250-300 тысяч (для сложных). При- 


мер расчетной сетки возле балки трапециевидного сечения с перекрытием и ограждениями приведен на рис 1. 


Рис. 1. Пример расчетной сетки (фрагмент) 


При решении уравнений Навье — Стокса связь скорость — давление реализовывалась с помощью ал- 
горитма 5 1МРГЕ. Конвективные и вязкостные члены уравнений движения жидкости и переноса турбулентных 
параметров аппроксимировались схемами второго порядка точности. 

Численное интегрирование осуществлялось неявной по времени схемой 2-го порядка точности. Шаг 
интегрирования ДЕ составлял (0,02... 0,04) Н/У» (У, — скорость набегающего потока), т. е. при сходе вихрей с 
безразмерной частотой ДН /У., = 0,1 был приблизительно в 250—300 раз меньше периода 1/{, и это обеспечи- 
вало приемлемое разрешение нестационарных параметров потока. Установившаяся вихревая дорожка обычно 
формировалась к моменту (60...120)Н/ У. Таким образом, общее число шагов интегрирования составляло в 
среднем 6000-10000. Для сбора нестационарной статистики использовался отрезок времени, не меньший 5 пе- 
риодов. 

Пример качественного сопоставления расчетной и опытной картин течения вблизи пролетного строе- 


ния приведен на рис. 2. 
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6) 
Рис. 2. Пример картины обтекания пролетного строения: эксперимент (а), расчет (6) 


Результаты исследования. Подробные сведения по АДХ типовых поперечных сечений можно найти 
в [11]. На рис. 3-14 представлено сравнение расчетных и опытных данных по коэффициентам осредненных 
аэродинамических сил (лобового сопротивления, подъемной силы) и момента для некоторых характерных се- 
чений. 





Рис. 3. Коэффициенты осредненных аэродинамических сил. 
Здесь: Сха — лобовое сопротивление; Са — подъемная сила; Си — момент; 


Ви Н — продольный и поперечный размеры сечения (без учета ограждений); а — угол атаки 


Механика 


Набегающий гладкий поток характеризуется интенсивностью 0,5 %, турбулентный — 8 %. Расчетные 
данные представлены сплошными линиями. 
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Рис. 4. Узкая одиночная прямоугольная балка (соотношение В/Н 0,75): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 5. Узкая одиночная прямоугольная балка (соотношение В/Н 1,24): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 


Механика 


67 


Вестник Донского государственного технического университета. 2016. Т. 18, № 4. С. 362-378. 155М 1992-5980 е155М 1992-6006 
Уезник о} оп Экие Тесйтса! Отуъегзйу. 2018. Уч. 16, по. 4, рр. 362-378. 1985М 1992-5980 е155М 1992-6006 

















В/Н = 2,3 
а) 





- Р-— 55 а8-0-0-542-0 9 0 9 


ы\я ооо || |] 
ОО 





Бир://уезбК.Аопза.га 


Рис. 6. Узкая одиночная балка с перекрытием (соотношение В/Н 2,3): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 7. Узкая одиночная балка с перекрытием и ограждениями (соотношение В/Н 2,3): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 369 
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Рис. 8. Широкая одиночная прямоугольная балка (соотношение В/Н 1,85): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 9. Широкая одиночная балка с перекрытием (соотношение В/Н 3,9): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 37] 
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Рис. 10. Трапециевидная балка (соотношение В/Н 3,09): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 11. Трапециевидная балка с перекрытием (соотношение В/Н 5,1): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 12. Двухбалочное строение (соотношение В/Н 2,75): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 13. Двухбалочное строениес перекрытием (соотношение В/Н 3,85): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Рис. 14. Многобалочное строение с перекрытием (соотношение В/Н 9): 
форма поперечного сечения (а); гладкий поток (6); турбулентный поток (в) 
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Обсуждение и заключения. Анализируя результаты, можно отметить, что за некоторыми исключени- 
ями с увеличением относительной ширины В/Н сечения однобалочного строения точность расчетного предска- 
зания его АДХ несколько возрастает. Как правило, наилучшее согласование отмечается по лобовому сопротив- 
лению сечения. Для большинства конфигураций расчетные данные несколько выше коэффициента лобового 
сопротивления, полученного экспериментально. Поясним, что для рассматриваемых плохообтекаемых тел пре- 
обладающий вклад в лобовое сопротивление сечения вносит сопротивление формы (давления), в основном 
определяемое разностью сил давления на наветренной и заветренной сторонах сечения. Принятый расчетный 
подход огрубляет динамику вихревых структур в области за телом, что приводит к заниженному восстановле- 
нию давления в указанной области. 

Величина подъемной силы более чувствительна к наличию, протяженности и типу (открыта/замкнутая) 
отрывных зон. Особенно это касается оснащенного плитой пролетного строения — в этом случае возможно 
повторное присоединение потока к верхней стороне плиты с образованием замкнутой зоны отрыва (примерно 
при В/Н >5). Поэтому по сравнению с лобовым сопротивлением расчетное определение подъемной силы ме- 
нее точно, особенно для пролетных строений с плитой перекрытия. 

Наиболее проблемным для большинства конфигураций является воспроизведение влияния угла атаки 
на аэродинамический момент сечения. 

Если при поперечном обтекании многобалочных пролетных строений между балками возникает аэро- 
динамическая интерференция [12], точность предсказания АДХ падает с увеличением числа балок (относитель- 
ной габаритной ширины сечения). В этом случае вместо ОКА\$-подхода целесообразно использовать более 
точные вихреразрешающие методы РЕЗ и [Е$. 
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